The Boson Peak and its Relation with Acoustic Attenuation in Glasses 
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Experimental results on the density of states and on the acoustic modes of glasses in the THz 
region are compared to the predictions of two categories of models. A recent one, solely based on 
an elastic instability, does not account for most observations. Good agreement without adjustable 
parameters is obtained with models including the existence of non-acoustic vibrational modes at 
THz frequency, providing in many cases a comprehensive picture for a range of glass anomalies. 
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The boson peak is an excess in the vibrational density 
of states (VDOS), g(uj), observed in many glasses at fre- 
quencies lv/2tt of the order of one THz. It appears as a 
hump in g{u>)/oj 2 vs. to, above the acoustic Debye level 
gv(u})/ui 2 . It is typically located at fiep ~ 0.1 x wd, 
where wd is the Debye frequency. This excess produces 
the well-known specific heat anomaly of glasses at tem- 
peratures T ~ M^bp/5£;b ~ 10 K [1]. It is generally 
agreed that the boson peak is a vibrational signature of 
the disordered structure of glasses beyond the nanometer 
scale. Its correct understanding is thus of considerable 
importance. Two main categories of dynamical boson- 
peak models currently exist. The first, which we call 
harmonic random matrix (HRM) models, is based on the 
concept that purely harmonic elastic disorder produces 
an excess of low frequency modes. The alternate picture 
is that there exist in glasses additional -non-acoustic- 
quasi-local vibrations (QLVs) at low frequencies |2|, |3(. 
The purpose of the present Letter is to compare HRM 
and QLV models to actual experimental results. We find 
that QLVs, using only independently determined param- 
eters, mostly provide a much better agreement between 
theory and experiment than HRMs can do. 

The HRM models postulate randomly fluctuating 
spring constants Ky, as e.g. in [J, la, la, 111 III- These 
assume a distribution p(K), generally bounded between 
two values K m i a < K max . A truncated gaussian distri- 
bution centered at K > but extending to K nlin < is 
used in [4| . The system becomes mechanically unstable if 
— K; min is above a threshold. Just before this instability, a 
low frequency excess in the VDOS appears asprecursor, 
which is interpreted as the boson peak [4j . In [5(] , a square 
distribution p(K) is used instead. A negative K m [ n is not 
necessary if pCK) oc 1/K over the interval {K m [ n , K max } 
as shown in [a]. Such a distribution can be rationalized 
by free- volume considerations. A similar distribution, in- 
cluding some negative force constants, was used in [7|] and 
handled with so-called Euclidean Random Matrices. The 
shape and position of the predicted boson peaks are not 
universal as they depend on the selected IVdistribution. 
While earlier treatments were based on microscopic har- 



monic models, more recently macroscopic tensorial elas- 
tic approaches were developed [a l2( • In [9( , the gaussian 
disorder affects only the shear modulus which tends to 
zero at critical disorder. This now allows performing re- 
alistic comparisons with experimental results, and it will 
be used below. In all these models, the boson peak is 
produced by shifting the maximum in g(oj)/uj 2 to low ui 
owing to a softening of the elastic constant distribution. 

On the other hand, inelastic neutron-scattering mea- 
surements of the dynamic structure factor of silica 
showed unambiguously that its boson-peak modes are 
not sound waves [id ], providing an early justification for 
the QLV hypothesis. That the boson peak of this net- 
work glass relates to librations of structural units has 
since been confirmed by hyper- Raman scattering 111 ], 
and similarly for boron oxide [l2j ■ Additional manifesta- 
tions of excess excitations are the thermal and acoustic 
anomalies observed below liquid-He temperatures that 
are described by two-level systems [1.3]. A theory en- 
compassing these and QLVs is the soft potential model 
[3,[lJ|. The latter predicts an onset in the excess VDOS, 
with gy/u> 2 oc ui 2 14 1. evolving into a 



gv = g - 3d, with gv/uj- oc lu- 
boson peak when the total number of QLVs is appro- 
priately limited [151 ] . This is now further understood 
on the basis of a physical model that considers the me- 
chanical instability of an initially uniform distribution of 
quasi-harmonic QLVs of density go(io), interacting via 
their common strain field [la, [l7| . The authors assume 
that go(oj) has a high cut-off, at ujq ~ u>d- The inter- 
action being weak, it only destabilizes those oscillators 
of relatively low frequencies, below uj c <C ojq [18]. The 
destabilized oscillators, restricted by their anharmonic- 
ity, form a renormalized density g c (oj) oc a; up to u> c [lqj . 
This result is independent from the initial go(uj) as well 
as from the size of the anharmonicity. A last but cru- 
cial step is that, owing to anharmonicity, the displaced 
oscillators produce random static forces affecting other 
nearby oscillators. This creates a soft gap, <?v{w)_oc w 4 , 
up to uj\y < uj c . The reader is referred to [la, [l7| for 
details. This model thus leads to a peak in gy(uj)/oj 2 , 
with an onset in u> 2 up to ~ Wb> and a decay in ss lo^ 1 
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FIG. 1: The measured boson peaks of ii-Si02 [21| and glycerol 
[23 ] compared to the Debye levels (dashed horizontal lines), 
to the best prediction of the HRM model [S| (dotted lines), 
and to the onset calculated using the soft-potential model, 
Eq. (la) (thin solid lines) and Eqs. (la+b) (thicker solid 
lines), (color online) 



between lo^ and lo c . While the onset power is univer- 
sal, the decay above Wb depends on the separation be- 
tween Wb and u) c . For small separation, the decay is be- 
tween w _1 and uj~ 2 [17] . in good agreement with obser- 
vation 19]. Contrary to HRMs, the QLV boson peak is 
not a shifted down end of acoustic branches. The high 
frequency VDOS is but slightly modified by the succes- 
sive interactions. Complemented with the soft-potential 
model, the QLV model makes specific predictions con- 
cerning high-frequency acoustic modes and their Ioffc- 
Regel limit beyond which they cease existing as plane 
171 [20| . Those will be checked below. 



waves 

First consider the VDOS in the boson peak region. 
High quality data are available on many glasses, in par- 
ticular from neutron scattering. For example, the exper- 
imental results, normalized to unity and divided by E 2 , 
where E = hui, are shown in Fig. 1 for the strong network 
glass silica, v-SiC>2 [21], and for the intermediate molec- 
ular one glycerol [22|. The corresponding Debye levels, 
g-£,(E)/E 2 , are shown by dashed lines. These are calcu- 
lated from the known velocities of longitudinal (LA) and 
transverse (TA) acoustic modes, «l and «t respectively. 
The boson peak of V-S1O2, at -Ebp = MIbp — 4 meV, 
exhibits the well known large relative excess gv / 9t> ~ 5. 
In glycerol, this ratio is much smaller, ~ 1.2. The dotted 
lines in Fig. 1 illustrate for comparison the predictions 
of the HRM model, Eq. (2) of [9]. These are determined 
itcratively using the experimental values of the sound ve- 
locities and of fiBp/^D to extract the appropriate bare 
velocities and separation parameter e of [9J] . It is obvious 
from Fig. 1 that the elastic instability model falls short 
of reproducing the observed peak strength. This already 
suggests that additional modes should be involved. For 
further comparison, the prediction of soft potentials is 
calculated below £bp , in the region where the growth of 
gv/E 2 is in E 2 , using [HOI] 
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the eigenvalue zero, calculated per atom, and W is 
the crossover energy between vibrational and tunneling 
states. These parameters are obtained from thermody- 
namical data around liquid- helium temperature (T) , and 
are available, both for silica [23| and for glycerol [24|. Us- 
ing (la), and adding the Debye contribution gu(E)/E 2 , 
one obtains the thin lines shown in Fig. 1. It is well 
known that there exist, in addition to the VDOS, quasi- 
elastic contributions to the scattering signal arising from 
relaxations in double-well potentials. This is described 
with the same soft potential parameters [l9|], 
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where k^T is the thermal energy. Adding g rc \ to gy, a 
quantitative agreement is obtained in the onset region be- 
low -Ebp , shown by thick lines in Fig. 1. This is not a triv- 
ial result, as the parameters P s and W entering Eqs.(l) 
are derived from independent very low T measurements. 
It emphasizes the self-consistency of QLVs and two-level 
systems embodied in the soft potential model. It is an 
additional indication that the excess indeed arises from 
QLVs. The solid lines in Fig. 1 do not account for the 
saturation in the growth of g(E)/E 2 near Ebp- This 
could be easily included, as shown e.g. in [19j]. According 
to 15|, [l7 1 , it requires at least one additional parameter 
which is beyond the purpose of the present discussion. 

A second important information is the linewidth of 
acoustic modes of very high frequency £l/2n, which near 
the THz range should increase as T ex O 4 , both in the 
HRM and QLV models. In HRMs, the acoustic width 
arises from elastic Rayleigh scattering by disorder. With 
soft potentials, it is the resonant absorption of sound by 
excess modes that gives the leading contribution [23| , 
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Here P s is the density of additional modes around 



Here, F is the full width of LA modes relating to the 
energy mean free path £ — Vi,/T, E is the phonon en- 
ergy hil, and Cl is the tunneling strength for LA modes 
that can be obtained from the height of the acoustic at- 
tenuation plateau near liquid-He T. This strong growth, 
r oc E" 4 , leads to a rapid end of acoustic plane waves 
as E increases. This occurs at a Ioffe-Regel crossover 
Em = fifiiR where £ decreases down to half the wave- 
length, i.e. for T = Tib, = ^ir/V [2^]. This crossover 
happens to fall at frequencies and wavevectors well above 
those reached in optical or UV Brillouin scattering, but 
most often near the lowest values reachable with inelastic 
x-ray scattering. The energy Em can be determined ex- 
perimentally provided T(E) is available over an adequate 
range [25(. It can also be calculated from Eq. (2) and the 
known soft-potential parameters, Em = 2(tt 2 Cl)~ 1 / 3 W. 
In the few cases where both can be done, their agreement 
is remarkable 12611. 
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FIG. 2: The full widths ftF of x-ray Brillouin spectra mea- 
sured on nine glasses in function of the mode energy E, both 
derived from adjustments with damped harmonic oscillator 
lineshapes [H, [27|, H, [H, El El El El El- The dotted 
curves are calculated following |9| as explained in the text, 
leading to the parameter e which is indicated. The lines of 
slope 4 are either calculated using Eq. (2) (solid lines), or 
traced through the experimental points (dashed lines), ter- 
minating then at the Em values obtained from the data as 
explained in [251 ] . Arrows indicate the approximate boson- 
peak positions taken from [35| where the experimental J7ir 
values are also discussed, (color online) 



The current situation for nine glasses is summarized in 
Fig;_2._Thedots show the x-ray data [H, [27|, [H, [H, 0, 
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32|, [33j, [34| . These are full widths of damped harmonic 



oscillator lineshapes adjusted to x-ray Brillouin spectra 
after convolution with the instrumental response. The 
HRM predictions (dotted lines) are calculated using Eq. 
(4) of [9j . Again the known values of vl,vt, and flsp /wd 
are used to extract the separation parameter e. It is clear 
that in all but one case the best HRM predictions fall well 
below the observed widths. Interestingly, the single ex- 
ception is a metallic glass [34f for which a real Ioffe-Regel 
crossover was not observed [361 ] and might not exist. For 
comparison, the QLV predictions below the crossover are 
drawn as solid lines for the three cases for which reli- 
able values of both W and Cl are available 37|. Above 
the crossover, Eq. (2) ceases to apply, and the widths 
cease increasing in E A . This presentation shows that the 
experimental widths smoothly prolongate the QLV pre- 
diction, which again supports the quantitative validity 
of this model on the basis of independently determined 
parameters. However, it also illustrates that the x-ray 
data unfortunately starts at somewhat too high energy 
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FIG. 3: Two inelastic spectra obtained on lithium diborate 
[25| at wavevectors Q near and above the Ioffe-Regel crossover 
at Qir ~ 2.1 nm~ . The dotted curves are the best predic- 
tions of the HRM model [S| . The relative amplitudes are sig- 
nificant. The instrumental full widths at half maximum for 
the Brillouin peaks, including the effect of the finite collection 
angle, are indicated by horizontal bars, (color online) 



in these three glasses to properly investigate the onset 
region ff oc £ 4 . On the opposite, this onset falls at suf- 
ficiently high E to be well observed in lithium diborate 
LI20-2B203 12 5| and in densified silica glass of density 2.6 
g/cm 3 , <i-Si02 [33j. In these cases, as well as for the four 
other glasses in Fig. 2, either one or both soft-potential 
parameters are unknown. For glycerol and polybutadi- 
ene W is known, but for C one only has an average value 
derived from low T thermal conductivity. It turns out 
that Cl and Cr (for TA modes) can differ by factors as 
large as ~ 3. Thus, the solid line in Fig. 2 cannot be 
drawn for these two glasses. Returning to the HRM pre- 
dictions, it is also remarkable that the slope of the dotted 
lines tapers off gently in all cases, and that no single line 
crosses the value T — Q/tt within the ranges shown. This 
means that the HRM model does not predict Ioffe-Regel 
crossovers at places where these are in fact observed [35| . 

It is also of interest to compare spectral shapes with 
specific HRM predictions of [9] . This is illustrated in Fig. 
3 for spectra observed on Li20-2B203 [25|. In drawing 
Fig. 3, the HRM lineshapes have been convoluted with 
the instrumental response function, and their height was 
adjusted to the corresponding experimental peak ampli- 
tude at low Q < Qib,. The instrumental response being 
wider than the spectral broadening, the latter is hard to 
judge just inspecting Fig. 3a. Only the numerical analy- 
sis of the data leads to the widths given in the correspond- 
ing panel in Fig. 2. However, as O increases beyond fiip,,, 
the HRM predictions, including the peak position, start 
departing strongly from the observed spectra which also 
become very wide. This presentation clearly shows that 
the best HRM prediction cannot provide an explanation 
for the observed spectra. In the QLV model, one expects 
that in the crossover region the acoustic modes hybridize 
with the boson-peak modes. In such a strong coupling 
case, perturbation approaches become inadequate. For 
this reason, a reliable expression for the Brillouin line- 
shapes in the QLV model is not yet available. 



In conclusion, it clearly appears that acoustic modes 
alone are not able to account for observations on most 
glasses in the THz range. The acoustic modes are merely 
a subset of all the low frequency vibrations in disordered 
systems. QLVs predicted by the soft potential model take 
these into account using parameters that arc determined 
fully independently by two-level-system measurements 
[13,1141]. It is truly remarkable that the latter, performed 
at low frequencies and sub-liquid-He temperatures, agree 
so well with hypersonic measurements at elevated tem- 
peratures. For network glasses, there exists separate ev- 
idence on the non-acoustic nature of the boson-peak ex- 
cess which relates to rigid librations of structural units 
IfJ, [n|, LL2J • For molecular glasses, the disagreement be- 
tween the observed widths and the HRM predictions in 
Fig. 2 is also striking. It is intuitive that for these the 
lowest frequency non-acoustic modes could be rigid libra- 
tions as well, though there exists so far little experimental 
evidence for this. Contrary to crystals in which modes are 
orthogonal, in disordered systems the QLVs can linearly 
couple to strains and to planar acoustic modes which are 
not true eigenmodes. For this reason the acoustic modes 
seen in Brillouin scattering are strongly affected as their 
frequency nears the boson peak. The observed relation 
between E\r and Ebp [20, I25l | 35l| is also well explained 
by this coupling, as shown in [171 ]. 
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These are obtained from the tunneling value C L 

iTUN , 



using 



the relation C L ~ C£' UN /1.1 [15(. For w-Si0 2 and Se, 



C T TUN is from J.F. Berret and M. Meissner, Z 



j is num j.r. oeneL anu ivi. ivieissiier, Li. Phys. B 

70, 65 (1988), for B2O3 it is from M. Devaud et al, Solid 
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